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In order to investigate the metallic state in high-mobility Si-MOS structures, we have further 
developed and precised the dipole trap model which was originally proposed by B. L. Altshuler and 
D. L. Maslov [Phys. Rev. Lett. 82, 145 (1999)]. Our additional numerical treatment enables us to 
drop several approximations and to introduce a limited spatial depth of the trap states inside the 
oxide as well as to include a distribution of trap energies. It turns out that a pronounced metallic 
state can be caused by such trap states at appropriate energies whose behavior is in good agreement 
with experimental observations. 
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I. INTRODUCTION 

The discovery of the metal- insulator transition (MIT) 
in two-dimensional (2D) electron systems in 1994^''^ has 
attracted large attention, as it was in apparent con- 
tradiction to the scaling theory of localization^i^ which 
states that in the limit of zero temperature, a metal- 
lic state should exist only in three dimensional systems, 
whereas in two dimensions disorder should always be 
strong enough to induce an insulating state. The MIT 
in high-mobility n-type silicon inversion layers shows a 
strong decrease of resistivity p towards low temperature 
T for high electron densities, manifesting the metallic re- 
gion, whereas a strong exponential increase of p demon- 
strated the insulting regime for low densities. A similar 
but weaker behavior was observed in many other semi- 
conductor systems at low densities and low temperatures 
(e. g. p-GaAs^ n-GaAs,^ SiGe,I AlAsS) 

Several models were suggested in order to explain 
the unexpected finding of metallic behavior in 2D. 
The most important ones are i) temperature-dependent 
screening)2i"— ii) quantum corrections in the diffusive 
regimefi^^— and iii) quantum corrections in the ballistic 
regimeiiSii^ Numerous argumentations for the different 
models are given in literaturej^ir— but a clear decision 
for one of them could not been drawn yet. 

As an alternative, Altshuler and Maslov (AM) intro- 
duced the dipol scattering scenario for Si-MOS struc- 
tures in which charged trap states in the oxide layer form 
dipols together with the image charge of the screening 
2D electrons3^ The interplay between the gate voltage 
dependent energetic position of the trap states and the 
height of the chemical potential may lead as well to a 
metal-insulator transition in that system. It should not 
be assumed that the dipol scattering effect is active alone, 
as the temperature dependence of screening and quantum 
corrections will surely contribute at low temperatures, 
but the charging of trap states might be the generator 
of the particularly large effect in Si-MOS structures. It 



is known that the misfit at the silicon/silicon-oxide in- 
terface produces charged defect states in the thermally 
grown oxide layer Arguments on the importance of 
trap states in Si-MOS structures were also given by Klap- 
wijk and Das Sarma.^° 

AM could show within their analytical calculations 
that a trap level at energy Ex which is either filled or 
empty, depending on its position relative to the Fermi 
energy Ep, can lead to a critical behavior in electron 
scattering if Et and Ep are degenerate. This dipole trap 
model is able to explain the main properties of the metal- 
insulator transition in gated Si-MOS structures. 

For the analytical calculations AM made a number of 
assumptions. These are: 

al) the trap states possess a (5-like distribution in en- 
ergy (i. e. have all the same energy), 

a2) the spatial density distribution in the oxide is ho- 
mogeneous, 

a3) the states occupied with electrons behave neutral 
and cause no scattering of 2D electrons whereas the un- 
occupied states are positively charged and lead to scat- 
tering (AM work in the hole trap picture, but we describe 
occupation in terms of electrons), 

a4) a charged trap state is screened by the 2D electrons 
so that the resulting electrostatic potential can be de- 
scribed by the trap charge and an apparent mirror charge 
with opposite sign on the other side of the interface, 

a5) the scattering efficiency of the 2D electrons is de- 
scribed by a dipole field of the trap charge and its mirror 
charge, 

a6) a parabolic saddle-point approximation for the ef- 
fective potential between the Si/Si02 interface and the 
metallic gate was used in order to perform analytical cal- 
culations, 

a7) the energy of the trap state Et is fixed relative to 
the quantization energy Eq of the 2D ground state inside 
the inversion potential, and 

a8) the chemical potential ji in the 2D layer has (A) 
either the same temperature dependence as in the bulk 
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substrate or (B) as in a 2D electron system with constant 
electron density. 

In this work, we precise and develop the AM trap 
model further in order to better understand the influ- 
ence of charged traps on the metallic state in Si-MOS 
structures. We present detailed numerical calculations of 
the temperature and density dependent resistivity p due 
to electronic scattering in the dipole trap model. Due to 
the numerical treatment we were able to drop the approx- 
imations a6), a7), and a8) of the analytic AM model. In 
addition, we have further extended our calculations for 
the more realistic case with energetic broadening and spa- 
tial distribution profile of the defect states, i. e. dropping 
also approximations al) and a2). 

As a result of our calculations, we find good agreement 
with the calculations of AM. There are mainly deviations 
in the overall behavior of the resistivity at low electron 
densities and at high temperatures. In order to under- 
stand the approximations of AM, we have also recalcu- 
lated the analytical model and were able to formulate it 
in a simplified way. Several mathematical terms are rear- 
ranged so that the scattering efficiency is expressed in the 
same form as in the usual Drude formulation. In addi- 
tion, their result contains two integrals, which we could 
replace by Fermi-Dirac integrals. Thus the known ap- 
proximations for the Fermi-Dirac integrals lead to simple 
equations for the resistivity p at low temperatures. 

The paper is organized as follows. In Sec. |TT] the an- 
alytical formulation of the dipole trap model is given in 
detail and we show that within the saddle-point approxi- 
mation the result can be written in terms of Fermi-Dirac 
integrals. Section IIIII treats the chemical potential and 
Sec. IIVI the analytic approximations for low tempera- 
tures. From the numerical integration, first results are 
given in Sec. |V] whereas in Sec. IVII the calculations are 
extended for the case that the conduction band is the ref- 
erence energy for the trap energy and not the electronic 
ground state Eq. In Sec. IVIII the model and the calcula- 
tion are extended for a spatial density distribution of the 
trap states and in Sec. IVIIll energetic broadening of trap 
states is taken into account. Conclusions are drawn in 
Sec. IIXI In two appendices the behavior of the chemical 
potential and the ground state energy of the inversion 
layer are described in detail. Please note that we use SI 
units throughout this work. 
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Figure 1. Schematic representation of the trap states together 
with the 2D electron system in the Si inversion layer. For 
Et > p, {T) the trap state is positively charged and scatters 
electrons in the 2D layer whereas for Et < p (T) the trap is 
neutral. 



denotes ^(T = 0). 

A potential gradient due to an applied gate voltage 
Vg causes a decrease of the trap energy Et = Ets — 
eVinsZ/ D, with the unscreened trap energy Ets at the 
oxide semiconductor (OS) interface, the voltage drop 
across the oxide (insulator) Vins, the distance from the 
OS interface Z, and the thickness of the oxide D. In 
their corresponding equation AM use the symbol Vg in- 
stead of V^ins- • But the total gate voltage Vg is equal 
to the voltage drop across the oxide Vins plus the voltage 
drop across the depletion layer Vdepi- Later in their paper 
AM write that the threshold voltage is incorporated into 
Vg, nevertheless they use practically an equation equiva- 
lent to 



Vr. 



en^D/einseo, 
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where e^s is the dielectric constant of the oxide and eg 
is the electric field constant. But if the threshold volt- 
age is incorporated into Vg, the latter cannot be used to 
calculate the slope of the electrical potential within the 
oxide, as a part of Vg falls off between the OS interface 
and the bulk layer. As the charges within the depletion 
layer (2D charge density — endcpi) also contribute to the 
gradient of the potential, we use the equation 



Vins ^ e{ns+ Udcpi) D/eins£o 



(2) 



together with 



II. TRAPMODEL 

In the AM model it is assumed that a large number 
of hole trap states exists in the oxide at a certain trap 
energy Et- If the trap energy lies above the chemical 
potential p, the trap is empty (in the electron picture, 
or has captured a hole in the equivalent description) and 
is positively charged, whereas if Et lies below p it is 
filled with an electron and thus is neutral, see Fig. [T] 
Please note that we use in this work the terminology 
p{T) for the chemical potential, the Fermi energy Ep 
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According to AM another term has to be added to the 
trap energy Et which describes the interaction with the 
two dimensional electron gas (2DEG) because of the im- 
age force. In Ref. [26l this term is given as — e^/STreinseo-^ 
(translated to the SI unit system). We think there 
should be a factor 16 in the denominator, see for instance 
Ref. [23 or Ref. HI- This factor does not change the re- 
sults of the trap model qualitatively therefore we write 
-e'^/CneinseoZ. 

AM define an energy 
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tering cross section at of this dipole they found classically 
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where Z is the distance between the trap and the oxide 
semiconductor interface, e is the energy of the scattered 
electron relative to the ground state energy of the inver- 
sion layer Eq: 



e — E — Eq, 
and e* is an effective dielectric constant 



(10) 



Figure 2. Oxide semiconductor interface. For simplicity we 
use different coordinate systems for the oxide (Z) and the 
semiconductor side (z) so that both, Z and z, axe positive 
on their sides. Trap energy Et which reaches its maximum 
a.t Z = Z,n, unscreened trap energy Ets at the OS inter- 
face, electrostatic energy Se, its components —eVinsZ/D and 
—edD/Z, and its maximum Sm, chemical potential fi, ground 
state energy of the inversion layer Eo and the corresponding 
wave function ip, conduction band edge Ec and its value at 
the interface Ecs- 



with C = 8 (AM) or C = 16 (our assumption). (We 
use in this work the notation of capital E for absolute 
energies and Greek e for energy differences.) Finally the 
trap energy can be written as 



Et (Z) = Ets + Se [Z) , 



ee{Z) 



D 



D 



(5) 
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where the subscript e in Eg stands for 'electrostatic'. The 
shape of Et [Z) is shown in Fig. [5] together with the 
relevant energy notation. 

Without the presence of a magnetic field the probabil- 
ity of a trap to be charged is given by 
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This formula is similar to the Fermi distribution function, 
differences are the minus sign in the exponent and the 
pref actor 1/2. A trap is charged when it is not occupied 
by an electron, thus the minus sign. If the trap is charged 
there are two possibilities for the spin orientation but 
only one if it is not charged, leading to the factor 1/2. 
As mentioned before p+ is determined by Et — /i (the 
vertical lines in Fig. |2l and the temperature. 

AM assume that a positive charged trap is screened by 
the electrons in the 2DEG and the trap forms together 
with that image charge a dipole. For the transport scat- 



(11) 



where esc is the dielectric constant of the semiconductor. 
We have recalculated Eq. (j5| and got the same result as 
AM. 

The Drude formula together with the Boltzmann equa- 
tion in relaxation time approximation yields the resis- 
tivity dp (Z) caused by the charged traps within the 
layer[Z, Z + dZ]. By integrating this contributions over 
the whole oxide [0, D] one gets 



V2rrr,c^e^/'' /n N+Z^^dZ 
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where rric is the conductivity mass of the free electrons 
within the inversion layer and e is an effective elec- 
tron energy as used by AM. In the corresponding AM 
equation the argument of the cosh function contains the 
Fermi energy, but should be replaced by the chemical 
potential. '^^ Furthermore Spo — Ep — Eq is the Fermi 
energy, fipo — fJ- ^ Eq is the chemical potential, each rel- 
ative to the ground state energy of the inversion layer, 
and is the 3D density of charged traps. With the 3D 
trap density Np [Z] it is given by 



N+{Z)^Nt{Z)p+ [Z). 



(14) 



Now we like to present Eq. in a form similar to 
the Drude formula 



P 



(15) 



and further find a term for the scattering rate 1 / (r) in 
the usual form that scattering rate is equal to scattering 
cross section times density of scattering centers times ve- 
locity of scattered particles, which leads to one of the 
basic equations used throughout this work 



mcCTt (e, Z) n+w (e) 



(16) 
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Here at (e, Z) is the transport scattering cross section 
from Eq. ^ for the effective electron energy e and an 
effective distance Z between the traps and the OS inter- 
face, = (Z) dZ is the 2D density of charged 
traps, and v {e) is the electron velocity which corresponds 
with £ = iricV^ /2. Comparing Eq. and ([TBj) we find 



N+ [Z) Z^/^Z _ N+ (Z) Z^'HZ 



^2/3 



^I^N+{Z)dZ 



(17) 



for this factor in at (e, Z) = Cae~^^^ Z'^^^ ■ 

In order to calculate the resistivity p the knowledge 
of is not necessary, cancels out with the denom- 
inator of Z^/^ within at (e, Z) , AM do not use it. As 
mentioned in the introduction we are interested in the 
metal-insulator transition depending on the electron den- 
sity rig, i.e. we want to know the temperature behavior 
of p as a function of rig- In this context is very useful 
in order to understand on the basis of Eq. that it 
contributes the main variations to the resistivity p [us , T) 
whereas at (e, Z^ and v (e) show only weak dependence 
on Us and T. The benefit of Eq. (|T6|) against (|T2|) is, that 
the physical meaning of the terms is immediately clear. 

The integral in the numerator and that in the denom- 
inator of Z'^/^ can be treated in quite the same way, so 
we define 
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N+ {Z) Z^dZ 
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(18) 



In the last step we followed AM and assumed that the 
trap density is constant within the oxide, respectively in 
the region where {Z) does not vanish. Now we can 
write 

at = c.£-i/3z2/3 = c.e--i/3^^ (19) 

"0 



p 



(20) 
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To be able to calculate the integral which corresponds 
with r22/3 AM expanded the electrostatic energy Et {Z) 
into a Taylor series about the point Z„i where it reaches 
its maximum Sm- This procedure is called saddle-point 
approximation. 
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Figure 3. Oxide-semiconductor interface. Trap energy Et (Z) 
(full line), its Taylor approximation (dashed line), and the 
resulting probabilities p+ (Z) in arbitrary units with p. as zero 
point. 



see Fig. |2] and |3l Now (fT8|) can be written as 



D 



Z^dZ 



I exp 



StsO = Ets - Eq^ const. 



(25) 
(26) 



AM assume that the energy Ets relative to the ground 
state energy Eq is constant, but we believe that rather 
the conduction band edge at the interface Ecs has to be 
used as reference energy, i. e. StsCs = Ets — Ecs — const. 
This issue will be further treated in section IVII 

The integrand is a peak around Z„i which drops off ex- 
ponentially on both sides. In order to come to the same 
result as AM, we further apply the following approxima- 
tions: (i) Because of the exponential decrease one can 
integrate from — oo to oo. (ii) The integrand is domi- 
nated by the denominator, so Z^ ~ can be set in the 
numerator. Now the integrand is symmetric around Zm 
and with help of the substitution Z — ^l^ifrp {Z — Z^f 
one gets 



z-^i^dz 



exp( ^"''7;;!^'^'" -In2 + Zl +1 



(27) 



This corresponds to the integral in equation (9c) in 
Ref. [2^ (Z corresponds to x^). We brought it into 
the form above as it corresponds now to a Fermi-Dirac 
integraP^ 



Z'^dZ 



r (A: + 1) Jo exp {Z - j]) + 1' 



(28) 
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where T is the Gamma function. A comparison yields 

x^^_,/,(^ln2+ ^^^^ J. (29) 

On the right hand side j appears only in the exponent of 
Z„i, so within the saddle-point approximation Eq. (|17p 
simplifies to 




^2/3 ^ ^ ^ ^2/3, (30) 



III. CHEMICAL POTENTIAL 

AM described two scenarios for the temperature be- 
havior of the chemical potential: (A) The chemical po- 
tential of the 2DEG and of the Si substrate coincide. (B) 
The 2DEG is disconnected form the substrate. For the 
case (A) they assumed that the temperature behavior in 
the 2DEG is the same as in the bulk. However, they did 
not take into account that the chemical potential in the 
2DEG is measured against the ground state energy Eq 
and in the bulk against the conduction or valence band 
edge, i. e. they assumed Eq and the band bending to be 
fixed. 

For (B) AM used an equation analogous to 



exp 



- 1 



(31) 



If only one subband of the inversion layer is occupied 
(quantum limit), the Fermi energy relative to its ground 
state energy is given by2S 



5s3v2D'Tld2D ' 



(32) 



where gs, gv2D, and md2D are the spin degeneracy, the 
valley degeneracy (for the 2DEG), and the density-of- 
states mass (2D) respectively. The two equations above 
can be derived from 



5s5v2DWd2D 

27rft2 



2D density of states 



1 



exp(f^)+l 



dE. (33) 



Fermi— Dirac distribution 



Our assumptions are shown schematically in Fig. |3] 
(semiconductor side of the OS interface) and are as fol- 
lows. In thermal equilibrium there is a single chemical 
potential ^ throughout the structure. For a certain tem- 
perature T, /i is determined in the bulk by the (residual) 
doping density (giving ^cb = A* ~ -E'cfc; for details see 
App. IB]) . In the inversion layer the position of /i relative 
to Eq (i. e. fjL^g) follows directly from the 2D density Ug. 



Figure 4. Band bending and notation on semiconductor side 
of the OS interface for thermal equilibrium between 2D layer 
and bulk. The electronic ground state energy of the inversion 
layer Eq, the conduction band edge Eq, its values at the in- 
terface Ecs and in the bulk Ecb, the valence band edge Ev, 
and the chemical potential /i are shown schematically, also 
the ground state energy relative to the conduction band edge 
£o, the chemical potential relative to conduction band edge in 
the bulk yLcb and relative to the ground state energy are 
noted. 



The band bending Ecb ~ Ecs adjusts so that 



E, 



Cb 



Ecs — £0 + fJ-Eo ~ A^Cb- 



(34) 



If the band bending is increased, the quantum well gets 
narrower and the ground state energy increases. Thus Eq 
itself is a function of the band bending, a self consistent 
calculation solves this problem (fix point iteration). 

Another possibility is that the 2D electron system is at 
low temperatures decoupled from the bulk substrate and 
no common chemical potential exists. We do not treat 
this case in this work, but expect a similar behavior. 

Fig. [S] shows the behavior of /i(T, rig) which is crucial 
for the understanding of the behavior of p{T,ns). The 
chemical potential relative to the ground state energy of 
the inversion layer fJ-EaiT,ns) increases for decreasing T 
and for increasing n^. 

The parameters, which were used here and for the fol- 
lowing calculations are collected in Table HI We assumed 
a {001} plane for the silicon surface. 



IV. APPROXIMATION FOR LOW 
TEMPERATURES 

For given values T and Us we calculate the resistivity 
p with help of the equations ([2T|) and (|29|) . Beside the 
explicit temperature dependence also the chemical poten- 
tial and the effective electron energy e are functions 
of T. 

As a first step we replace the integral in by a 
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Figure 5. Chemical potential ^ relative to the ground state 
energy Eq. a) Curves with constant n^, bottom up: ris = 
1,2,5,10,20 X 10" cm"", b) Curves for constant T, top 
down: T = 0, 50, 100, 150, 200, 250, 300 K. 



Table I. Values used for calculations. 

gs = 2 

ffv2D = 2 
5v3D = 6 

me = 9.1094 X 10"=*^ kg 
'7ld2D = 0.1905 me 
m^ = 0.9163 me 
mdo3D = 0.322 me 
mdhSD = 0.59 me 
esc = 11.9 
eins = 3.9 
e* = 7.9 
D = 200 nm 
C= 16 
ED = 0.4615 meV 
NA^'ix 10^^ cm~^ 

egosi = 1.17eV 
asi = 4.73 X lO-^'eVK"^ 
^Si = 636 K 



Fermi-Dirac integral. For the Fermi distribution function 
1 1 



exp (j^) + 1 exp 
the following identity holds, 

df{e) _ 



knT 



1 



(35) 



1 cosh-^fi-i^' 
2kBT 



de ikeT 
So Eq. can be written as 



-) 

K^FOJ 



de 



de 



(36) 



(37) 



An integration by parts yields 



e — epo 



1 -6 



£_F0 



-1/6 



f{e)de 
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6./0 



-1/6 



exp 

5 r 



Hyl-^-i/a 



knT 



-6 



n -6 



(38) 



For low temperatures we can take advantage of the 
behavior of fj,Eo (T) and e (T) for T — >■ 0. It can easily be 
shown that for constant Hs and therefore constant epo 
all derivatives vanish at this point. 



d"/XBo (T) 



rf"£(r) 
dT" 



T-fO 



for7i>l, 



for n > 1 . 



(39) 
(40) 



So /xfo (T') and e (T) are very flat functions at T — > 0, 
for ksT <^ epo they can be approximated by 



e{T)~eFo- 



(41) 
(42) 



See appendix [K\ for details. 

The Fermi-Dirac integral J-j (?/) can be approximated 
by22 



^1 iv) 



n 



r(j + 2) 



for 77 0, 



(43) 



which is the first term of an asymptotic series derived 
with help of a Sommerfeld expansion, '^^ and by 



Tj (77) ~ exp 77 for 77 <C 0. 



(44) 



For T — > the argument of the Fermi-Dirac integral in 
Eq. increases beyond any border. Which approxima- 
tion for the Fermi-Dirac iirtegral is applicable depends 
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on the sign of etsO + £m — £fo- (Here fiEa is replaced 
by Spo as by definition (j^Eq — t- £fo for T — > 0.) This 
is conform with AM's definition of the transition point, 
(^Tso + £m — £fo) /ksT = 0. Accordingly we define 

SmF > — > insulating, 

£mF < — > metallic, 

EmF = transition point 

and a critical density 

nsc = "-sle„^p=o • (45) 

Applying the appropriate approximations results in 

^ 2iVTZ^,+'/y ^ exp (l^) for£„^<0 
^-1/2 (ln2) A^T^i^+'/' for e„^F = 0. 

~0.891 

(46) 

Please note an interesting behavior. When setting SmF — 
in the first two equations for T > they converge nei- 
ther into each other nor into the third one. This apparent 
discrepancy can be understood, as |emF| gets smaller and 
smaller the maximum temperature where the approxima- 
tions for the Fermi-Dirac integral are just applicable also 
gets smaller and smaller and finally vanishes for SmP = 0. 
Indeed for EmF = and T — ^ the three cases yield the 
same result, i. e. flj — 0. 

V. COMPARISON OF ANALYTIC AND 
NUMERICAL RESULTS 



We like to explain the shape of p (T, Ug) starting from 
a value rig > 't^sc- The maximum of et (Z) lies below 
/i, therefore the (Z)-peak is very small and narrow. 
When ng is decreased /i drops off (see Fig. [S]) and the 
resistivity p rises fast because — {Et — p) /ksT is the ex- 
ponent in the denominator of . It rises the faster the 
smaller the temperature is. When p reaches the max- 
imum of the trap energy Et the p+-peak has a height 
of 2/3, therefore further decreasing of and p cannot 
increase the height of the peak appreciable (on a loga- 
rithmic scale, which covers some orders of magnitude), 
only the width. 

For increasing temperature T also the chemical poten- 
tial p measured against Eq decreases which results in 
increasing resistivity p. Additionally the p^ (Z)-peak is 
broadened because the Fermi distribution function de- 
clines over several ksT. 

The saddle-point approximation works best for small 
temperatures T and densities Us > rise, in this domain 
the (Z)-peak is very narrow and therefore the Taylor 
approximation of the trap energy is accurate within the 
peak. For large T and/or Ug < the deviations of the 
saddle-point approximation from the numerical calcula- 
tions is visible in Fig. [51 

It should be noted that according to the calculations, 
the resistivity p drops to arbitrarily low values in the 
metallic regime. This is caused by the fact that electron 
scattering is taken into account only by the trap states 
at a single trap energy. If this trap energy is below the 
Fermi energy, with decreasing T the number of charged 
scattering centers goes to zero. Only if other scattering 
effects, like residual impurities, surface roughness, accep- 
tor states in the depletion layer, etc. are included, the 
low-T resistivity would be limited. As will be seen later, 
an energetic broadening of the trap states will have a 
similar effect. 



In order to get rid of the restrictions which came from 
the saddle-point approximation we have performed nu- 
merical calculations of the integrals Q.j (Eq. pS]) ). In 
Fig. |S] we show the resistivity p depending on the tem- 
perature T and the 2D electron density in the inversion 
layer Ug calculated with the help of the saddle-point ap- 
proximation (full lines), the approximation for low tem- 
peratures (dashed lines), and the numerical integration 
(markers). We chose £tsO — 42.02 meV in order to get 
Use = 10"'^^cm~^^, for the 3D trap density we assumed 
Nt = 10^^ cm^'^. The value for Nt seems to be quite 
high, but only in a very narrow layer the traps will in- 
deed be charged (where the trap states are above the 
chemical potential) and contribute to scattering. A fur- 
ther limitation of the available trap states into a narrow 
region besides the interface will follow later in this work. 

The critical behavior of p (T) versus rig for the homo- 
geneous trap density Nt can be seen most clearly on the 
double logarithmic Figs. [6^) and c). We added [BId) as it 
is in this form directly comparable with Fig. lb) in the 
work of AM.2^ 



VI. CONDUCTION BAND AS REFERENCE 
ENERGY 

If the bands in the semiconductor and in the oxide are 
bent due to an applied gate voltage all energies move 
up and down with the bands. The energetic position of 
the trap states ETg should thus be fixed relative to the 
conduction band edge and not to the ground state energy 
Eq of the inversion layer as assumed by AM. 

From Eq. ^ we see that Et — P = ETg + Se {Z) — p 
determines the probability p^ of a trap to be charged. 
So if we measure ETg against Ecg we also have to know 
A*Cs = A* - Ecg- We find 

A*Cs ^ p- Eq+ Eq- Ecg = Pe„ + £q, (47) 

where Eq = Eq — Ecg is the ground state energy relative 
to the conduction band edge at the OS interface. 

The chemical potential pE„ follows from Eq. ([5T|l . but 
an accurate calculation of Eq is rather complex. For sim- 
plicity we follow the calculation method of Ando, Fowler 
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Figure 6. Behavior of the resistivity p depending on the temperature T and the 2D electron density in the inversion layer ris. 
Full lines: saddle-point approximation, point symbols: numerical integration, dashed lines: low temperature approximations. 
On the curves in a) - c) ris and in d) T is constant. Critical density rise = lO^^cm"^. a) logarithmic and b) linear display 
of p(r) with parameters: Us = 0.75, . . . , 1.25 x n^c in steps of 0.05 x rise in top down order for the individual curves, c) 
logarithmic view of p(r) with top down parameters: ris — 0.95, . . . , 1.05 x rise in steps of 0.01 x rise, d) p{ns) with parameters 
T = 0, 0.2, 0.5, 1, 2, 5, 10, 20 K in bottom up order. 



and Stern (AFS)^ and neglect the exchange interaction 
and correlation effects and use the Ritz variational prin- 
ciple. (In the mentioned article also more sophisticated 
methods for the calculation of Eq are given.) 

For convenience we introduce a new coordinate system 
z = ~Z, i. e. the z-axis is perpendicular to the OS inter- 
face, positive z-values correspond with the semiconduc- 
tor side. For the electrons in the inversion layer the bent 
conduction band of the semiconductor together with the 
step at the interface builds the quantum well. We use the 
Fang-Howard envelope wave function according to AFS^ 

/ ^,^ [A/¥^exp(-'f) forz>0 
ip{z,b) ^ ■{ \ ^ ^ ^ 2 ; - (48) 

1^0 forz<0. 

The parameter b is varied in order to make the total 
energy per electron minimal. For the potential several 
approximations are taken, see App. iBlfor details. 

In Fig. [7| the ground state energy Eq versus the elec- 
tron density rig is shown, Eq decreases with decreasing rig. 
Now we hold the difference between trap energy and con- 
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Figure 7. Ground state energy of the inversion layer eo versus 
the density of electrons in the inversion layer ris for three 
temperatures, top-down: T = 0, 200, 300 -ftT. 

duction band edge at the interface EtsCs — Ets — Ecs — 
Ets ^ + ^-Ea + ^0 constant (instead of EtsO ^-s before) . 



9 



When Us is decreased not only but also £o decreases, 
thus /i drops off faster against Ets and the transition 
is more abrupt. This can be seen in Fig. [5^) where re- 
sults for EtsCs — const, (full lines) and EtsO = const, 
(dashed lines) are compared. In both cases flj was cal- 
culated numerically. The critical curves — Usc coin- 
cide at T — because the values £tsO — 42.0 meV and 
EtsCs — 68.4 meV were chosen in order to get the same 
Use and for higher temperatures because {T) is nearly 
constant over a wide temperature range (Fig. [51:)). For 
the 3D trap density we assumed again Nt = 10^^ cm^"^. 

It has to be emphasized that we still use Eq. (HI). But 
the 2D density of positive charges in the depletion layer 
f^dcpi is a byproduct of the calculation of Eq. So it is no 
longer a problem to use Eq. ^ l^ns c>c -|- ndopi instead 
of ([T]) Vins ns- We will do so henceforward. As a result 
the slope of the line ETs — ^^ViasZ / D in Fig.[2]is increased, 
the maximum of Et (Z) falls off against /i, the (Z)- 
peak gets smaller and narrower. But it is still possible to 
let coincide the critical curves p{ns — Usc) by increasing 
the trap density Nt in order to compensate the narrower 
p+ (Z)-peak and by increasing StsCs iu order to get the 
same critical value Usc- Here it is also important that the 
2D charge carrier density of the depletion layer Udepi {T) 
has almost no dependence for low temperatures as 
can be seen in Fig. [SJi). The resistivity p{T,ns), with 
Nt = 2.98 X 10^*^ cm'^ and etsCs ^ 95.7 meV, for which 
we get the same Ugc as before, is shown in Fig. [BJa) . The 
above used depletion density ridcpi was calculated under 
the assumption of a background doping density of Nj^ — 
2 X 10^^ cm~'^, which is a typical value for high mobility 
Si-MOS samplesj^ Here the values for Nt and EtsCs 
were chosen in order to get the requested Usc for the 
given Nji. In reality the value EtsCs is determined by 
the chemical nature of the defect and thus the critical 
density may be different from sample to sample if the 
background doping is different. 

If now the slope of the energy Ets — eV-msZ/D 
is higher due to the inclusion of ridcpi, the variable 
max {Et {Z) — jj) which is crucial for the resistivity is 
less sensitive on Vi-^s aud Ug, therefore the transition is 
less abrupt. The Ug dependence of rtdcpi does not play a 
role because it hardly exists. 



VII. SPATIAL TRAP PROFILE 

At higher temperatures the large value of fc^T leads to 
charged trap states in regions where the chemical poten- 
tial is even somewhat below the chemical potential (i. e. 
the p+ (Z)-peak broadens) and therefore the resistivity is 
increased to unrealistic high values, see curves at higher 
temperature in Fig. [51 But an appreciable density of 
traps should exist only within the strained region of the 
oxide, and thus the broadening of the peak beyond the 
width of this region leads to an unrealistic description. 
We can resolve this problem by introducing a spatial trap 
density profile Nt{Z). 



If now the trap density Nt is a function of Z it has to 
remain inside the integral Q,j (compare Eq. I18|) 



I'D I'D 

I N+ {Z) Z'^dZ = / , 
Jo Jo 



Nt{Z)p+ [Z) Z'dZ. 



(49) 

For simplicity we use here an rectangular spatial trap 
profile from the OS interface to an arbitrary depth 2^max, 



Nt{Z)^1^^'' 



for < Z < Z„ 

forZ > Zmax, 



(50) 



where ut is the 2D trap density. Fig. [51 shows p{T,ns) 
for Zniax = 4 nm, the conduction band edge Ecs was 
used as reference energy, StsCs — 95.7 meV was chosen 
in order to get = lO^^cm"^. Where the 3D trap 
density Nt does not vanish its value is assumed to be 
2.98 X 10^* cm"'^ as before (see Fig. Hb)), resulting in 
TiT = Nt ■ Z„^i^x = 1.19 X 10^^ cm~^ of which again only 
a part is charged. 

As can be seen in Fig. [51 the behavior for low tem- 
peratures has hardly changed, but for high temperatures 
p (T, Tig) now saturates as a broadening of thep+(Z)-peak 
beyond .^max does not lead to a further increase in the 
number of charged scattering centers. This saturation of 
p{T,ng) is in fairly good agreement with experiments, 
where p for high T is limited as well. 



VIII. BROADENING OF THE TRAP ENERGY 

As the trap states will not all be identical and in ad- 
dition the stochastic position distribution will influence 
them mutually, their energetic position has to be broad- 
ened. 

We describe the broadening /S.Et with the help of a 
normalized distribution function g (^Ets, Ets, ^Et^ for 
the trap energy Ets which characterizes the trap, see 
Fig.[2| Now Ets has the meaning of a mean value. (Mean 
value should not be understood in a strict mathematical 
sense, e. g. for the Lorentz distribution the mean value 
does not exist, but in this case it is obvious to take the 
energy ETg where the distribution reaches its maximum.) 
Furthermore Ets is the value for a particular trap. The 
probability of Ets to lie within the interval [E, E + dE] 
is given by g {E, Ets, ^Et) dE. Therefore we replace the 
probability p+ of a trap to be charged by 



P+iZ) 



(^Ets , Et 



s,AEt 



-dEn 



(51) 

The denominator is that of p+ , only Ets is replaced by 
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Figure 8. Resistivity p (T, Us) for the potential across the oxide layer a) Vina oc Us with A'r = lO'^® cm~^ and stsCs — 68.4 meV 
and b) Vins oc + ridepi with Nt = 2.98 x 10^® cm~^ and etsCs — 95.7 meV. The parameters are chosen to give a critical 
density of Usc ~ 10^^ cm^^ for both cases. At the individual curves Us = 0.75, . . . , 1.25 x Usc in steps of 0.05 x n^c from top to 
bottom, full lines represent etsCs ~ const, (realistic case), dashed lines etso = const, (for comparison). In c) the ground state 
energy of the inversion layer eo {T,ns) is shown for same ris values as before, but now assigned bottom up. In d) the depletion 
density ridepi {T,ns) is shown, but the Us dependence vanishes within the line width. 




T(K) 



Figure 9. Resistivity p{T,ns) for rectangular spatial trap 
profile with Z^ax ~ 4 nm with curves of constant Us . The 



Ets- By introducing the dimensionless parameters 
AEt 



a = 



Ets — eVins 5 /U 



(52) 
(53) 



and a dimensionless distribution function h{ri) defined 
by 



7] = 



Ts ~ ^Ts 

AEt ' 



critical density is Usc ~ 10 cm 

has to be 95.7 meV. Densities for full lines are n 
0.75, . . . , 1.25 X rise in steps of 0.05 x rise, and for dashed 
lines Us = 0.96, . . . , 1.04 x rise in steps of 0.01 x rise (always 
top down), with Nt = ut/Z^^^ = 2.98 x 10^* cm"^. 



which means that etsCs g (ets, Ets, ^Et) = -i-^h (r?) , 

or full lines are n. = ^ ' AEt 



(54) 
(55) 
(56) 



the probability P+ can be written as 



P+(a 



/OO 
T 
-OO 2 



exp {—ar] — /3) + 1 



dr). (57) 
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As a rule this integral cannot be calculated analytically. 
An exception from this rule is the uniform distribution. 
If we define the width of the 'rectangle' as 2/S.Et we get 



and 



elsewhere 



p , 1 , l + 2cxp + 

F+ (a, «) = — In — r. 

+ 2a l + 2exp(/3-a) 



We also use the normal distribution 

1 / .2 



exp 



2 



(58) 



(59) 



(60) 



with the standard deviation as /S.Et and the Lorentz 
distribution (natural line broadening) 



1 1 



1 



(61) 



with the half full width at half maximum (FWHM) as 

Fig. [TOl shows the numerical results for the three differ- 
ent broadening distributions and different width l^E^ = 
0.02, 0.1 and 0.5 meV. Again we took the conduction 
band edge at the interface Ecs as reference energy, as- 
sumed that traps exist in the oxide only within 4 nm 
from the OS interface (constant trap density within this 
region), and used Eq. ^ instead of ((1]). We chose 
EtsCs — 95.7 meV in order to get n^c — 10^^ cm~^ and 
for the 3D trap density Nt — nT/Zma^ — 10^^ cm~^. 

As one would expect, the transition is less abrupt 
for higher AEt and does not vanish even for AEt = 
0.5 meV. If the broadening is indeed caused by potential 
fluctuations due to disorder or trap-trap interaction, the 
value of AEt could even be much larger than ksT and 
the transition might be smeared out even stronger. 

In the metallic regime the mean trap energy is below 
the chemical potential. As the normal and the Lorentz 
distribution have tails, there always remain some charged 
traps from the upper tail when otherwise all traps would 
be filled with electrons and therefore would be neutral. 
On a logarithmic resistivity scale the p (r)-behavior is 
changed drastically by the few additional charged traps. 

In the insulating regime the mean trap energy is above 
the chemical potential. Here a large part of the traps is 
charged and by contrast the few uncharged traps due to 
the lower tail of the energy distribution hardly play any 
role. 

By means of analytical considerations we got the fol- 
lowing estimations for the temperature below which 
the resistivity becomes almost constant. In general the 
resistivity p varies by some orders of magnitude, so as 
a criterion for being almost constant we took that re- 
gion where p changes finally by only a factor two down 
to zero temperature. The markers in Fig. [TO] represent 
these temperatures Tb. 



According to this definition, for the insulating behavior 
Us < Use for all three distributions we get 



maji Et (Z) - Ef 



(4 -In 2) fcs 
and for the metallic behavior > 



(62) 



AEi 



T — ) (4+ln2)/cB(BF-max£;T(2)) 
^ ~ \ Ep-mayj Et(Z) 



2(4+ln2)fc£ 



. normal distr. 
. Lorentz distr. 



(63) 

The uniform distribution has no tails so in the metallic 
regime there is no temperature range where p is almost 
constant. 



IX. CONCLUSIONS 

In this work we have performed numerical calcula- 
tions within the dipole trap model for Si-MOS structures. 
Originally this model was proposed by Altshuler and 
Maslov with several approximations, in order to be able 
to get analytical solutions. Due to our numerical treat- 
ment we could eliminate several approximations. We de- 
scribe the potential inside the insulator by its detailed 
spatial dependence instead of the saddle-point approxi- 
mation with the quadratic dependence around its max- 
imum, we fix the trap state energy relative to the con- 
duction band edge instead of relative to the electronic 
ground state inside the triangular potential well and we 
have taken into account the detailed change of the chem- 
ical potential in the two-dimensional electron layer with 
respect to the bulk material, which seems to be more 
realistic than the two cases in the original treatment. 

According to our calculations, the metallic regime at 
high electron densities n^, where the resistivity is de- 
creasing towards lower temperature, is strongly devel- 
oped . Also a critical density n^c can be identified with a 
characteristic temperature dependence different from the 
metallic and the 'insulating' region. For electron densi- 
ties Ug < Hsc, the resistivity curves satturate towards 
lower temperature and remain constant when the tem- 
perature approaches zero. They do not show an insult- 
ing behavior in the sense that p increases towards zero 
temperature. Such an increase can in principle be caused 
e. g. by a further decrease of the chemical potential p, with 
temperature, as in the work of Altshuler and Maslov for 
case A, where it was assumed that the temperature de- 
pendence of p ist the same for the two-dimensional elec- 
tron layer as it is in the Si-bulk material. In the dipol 
trap model a constant efficient screening is assumed in 
order to clarify the effects which are caused when the 
traps change their charge state. A realistic treatment 
of the temperatue dependence of the electronic screening 
could also cause an increase in p towards lower temper- 
ature in that it favors the formation of dipol trap states. 
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Figure 10. Resistivity p{T,ns) for different types of trap energy broadening functions and different width AEt- The picture 
columns show uniform distribution, normal distribution, Lorentz distribution from top to bottom. The rows show energy width 
of AEt = 0.02, 0.1, 0.5 meV from left to right. For the curves Us is constant. The critical density is rise = lO^^cm"^ for 
etsCs — 95.7meV. Full lines show rig = 0.75, . . . , 1.25 x n^c in steps of 0.05 x n^c , dashed lines: = 0.96, . . . , 1.04 x Usc in 
steps of 0.01 X rise (from top to bottom), for a 3D trap density of Nt = nT/.^max = 10^* cm~^. 



Furthermore, the quantum corrections in the weak and 
strong locaHzation regime are also neglected here, but 
would finally increase the resistance p at low T. 

In addition, we have further generalized the dipole 
trap model by dropping the assumptions that the trap 
states are homogeneously distributed inside the oxide 
layer and that the energy distribution is (5-like. A nar- 
row spatial distribution of the trap states near the oxide- 
semiconductor interface limits the number of charged 
states at high temperatures and thus gives an upper limit 
for the increase of the resistivity p as well. This leads 
to a good agreement with experimental observations at 
higher temperatures. The energetic broadening of the 
trap states on the other hand leads to a finite amount 
of unoccupied and thus charged states in cases where 
otherwise all states would lie below the chemical poten- 
tial /X and the number of charged trap states would go 



to zero for fc^T — > 0. Thus for high electron densities 
with metallic behavior the resistivity will not further de- 
crease towards lower temperature, but saturate at a finite 
values, as has been observed in experiments on Si-MOS 
structures as well. 

The effect of a magnetic field can be taken into account 
by the Zeeman splitting of the trap states with spin ±1/2. 
As shown by Althuler and Maslov, the energetic splitting 
can turn a metallic behavior into an insulating one. We 
did not include magnetic field effect in our calculations, 
but an according energetic shift of the trap states has to 
lead to the same effects in our refined model as well. 

We also like to mention that for low electron densities 
care has to be taken for the dipol scattering model. It is 
assumed that the electrons in the two-dimensional layer 
shield the potential of the charged trap states and thus 
form together a dipol field which is responsible for the 
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scattering. At very low electron densities this screening 
becomes weaker and the scattering will finally increase so 
that the resistivity should be higher in this regime. These 
effects have not been taken into account in the frame of 
the current work, as we like to present the basic effects 
due to charging of trap states. 

Althogether. our detailed numerical calculations 
within the dipol trap model show that a pronounced 
metallic state can be caused by trap states at an appro- 
priate energy level inside the oxide of Si-MOS structures. 
For the realistic assumptions of energetic broadening and 
narrow spatial distribution near the oxide-semiconductor 
interface, the behavior is in close agreement with exper- 
imental observations. 
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Appendix A: Temperature behavior of the chemical 
potential and of the effective electron energy for low 
temperatures 



First we show that 



= for 71 > 1 



ej?o=COnst., T^O 



(Al) 



holds for /i Eo {T) from Eq. ((3T|) . In our calculations we 
use the temperature T and the electron density in the 
inversion layer Us as independent variables. So it is al- 
lowed to set Ug = const, and therefore also Spo = const, 
while varying T, see Eq. ([5^ . 

For simplicity we introduce the auxiliary variable 



and the function 



M{x) 



£fo 



With the chain rule we find 



'FO 



So Eq. (|Aip is equivalent to 
d"M (x) 



for 71 > 1. 



da;" 

From Eq. (ISTI) we get 

M (x) — X In 



= f or 7T, > 1 . 



1 — exp — 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



The first two derivatives are 
dM 



dx 



In 



1 — exp 



1 



[exp(i)-l].' 
d^M _ 1 1 

^ "[cxp(i)-l] x-^ [cxp(i)-l]'x3- 



(A7) 
(A8) 



The second derivative contains only terms of the form 

1 



5a, 6 



[exp(i)-l]'^.^- 
Differentiating ga,b yields 



dx 



a5a,&+2 + a9a+i,b+2 - bga, 



6+1 1 



(A9) 



(AlO) 



therefore all higher derivatives 



also contain 



only terms gj_k with j G {1, 2, . . . , 7?.} and k g 
{ri + 1, n + 2, . . . , 2?! — 1}. Now we can multiply the nu- 
merator and the denominator of gj ^. with exp (— j'/x) and 
write 



T ( \ V exp(-|) 
lim Oa 1, (x) = lim ^ — — 



= lim 



lim 



exp I 



y->oo exp (y) ' 

with y — i/x. Applying the rule of L'Hospital k times 
we get 



lim 



fc! 



lim Oj,fe {x) - 
x->o y-^oo J**: exp (yj 



0. 



The limit for the first derivative also vanishes 
dM 



lim 

x-i-o dx 



In 1 — lim Qi 1 = 0, 



so Eq. (|A5|) holds and 

^iEo — EpQ for ksT ^ efo- 
With help of Eq. (jAl|) we show now that 



d"e(r) 
dT'' 



= for 7T, > 1 



£j?o=const., T-s-O 



(All) 



(A12) 



(A13) 



(AM) 



holds for Eq. 

For not too small 77 the Fermi-Dirac integral J^j (77) can 
be approximated bj*^ 



V 



J+1 



r(j + 2) 



(A15) 



For low temperatures jj^Eo approaches spo therefore 
it is positive, so from Eq. (j3H| we get 



( En 



ioikBT-^Epo, 



(A16) 
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and with the chain rule respectively Faa di Bruno's for- 
mula 



de 
dT 



de dfiEg 
d^Eo dT 



d^e 



. fc„! d^*-' 



k — ki 



1 dnEo 
1! dT 
■ ■ + fc„, 



£^0 
fci 



1 rf"Mgo 

7i! dT" 



where the sum runs over all integer numbers fci , 
which fulfill 



ki + 2fc2 



(A17) 



(A18) 



(A19) 



We do not have to find this numbers, we only need to 
know that for n > 1 at least one fc,- > 0, therefore each 



term in (|A18|) contains a factor 

T the sum vanishes and from Eq. (IA16P we get 



e^epo for ksT epo- 



^, so for n„ = const., 
;t 

(A20) 



Appendix B: Ground state energy of the inversion 
layer 

As mentioned before we calculate the ground state en- 
ergy of the inversion layer with help of the Ritz vari- 
ational principle using the Fang-Howard test envelope 
wave function^'' 



ip{z,b) 



exp I 



6z^ 
2 ; 



for z > 
forz < 



(Bl) 



and as an approximation for the potential 



U{z) 

Ud{z) 
Us{z) 

U^{z) 



Ud{z) + Us{z) + U,{z) iovz>0 
oo for z < 0, 



e rid, 



r-(l-2t^) iorz<Zd 



forz > Zrf, 



26esceo 
6e2 



6 - (^{bzf + Abz + cxp {-bz) 



327resceo Cs 



4z 



(B2) 

(B3) 

)" , 
(B4) 

(B5) 



The term Ud comes from the charged acceptors within the 
depletion layer with thickness z^. Us describes the inter- 
action with all other electrons in the inversion layer, and 
Ui the interaction with image charges. To write Us (z) 
in this form we have to assume that only the first sub- 
band is occupied, this is the so called quantum limit. 
The conduction band edge is built by Ud (z) -f Us (z), 
therefore the zero point of the energy scale was chosen 



to get Ud (0) -I- Us (0) = which means that the resulting 
ground state energy is measured against the conduction 
band edge at the interface Ecs as requested. 

The Hamiltonian is given hy H = T+U with the oper- 
ator for the kinetic energy T — ~ , where is the 

z-component of the effective mass of the semiconductor 
in the bulk. The ground state energy is the expectation 
value of the Hamiltonian, 



T 



U„ 



Us 



u. 



(B6) 



calculated with the value of b which makes the total en- 
ergy per electron 



T 



Ud 



1 



Us 



U, 



(B7) 



minimal. The factor 1/2 in the third term prevents from 
double counting the electron-electron interaction. With 
the density n* = ndcpi + ^^"^s introduced by AFS^S one 
gets 



-b^ +I3b + -lb" 



2 



4m, 



327resceo £1 



7 
5 



12e2 [Na - Nd) 




X n 



Na and No are the densities of the acceptors and donors 
respectively. The coefficients have been chosen in order 
to get a, /?, 7, i5 > and to get a most simple equation 



dr 

"g^ab + p- -fb-^ + S (b) b- 



0. 



(B14) 



As AFS used 



dcpl / 

£0 V 



1 



for < z < cx) instead 



of Eq. (jBSp they did not get the term proportional to 
exp {—bzd), under normal circumstances it is very small 
but formally it is necessary to see that b 0, e —00 
is not the global minimum. Furthermore they neglected 
/3 and S (b) b^^ and got 



6 = 



1/3 



12e n*m2 



1/3 



(B15) 
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We neglected only S {b)b ^, this leads to 



6: 
F ■ 
K 
B ■ 



F i^K 

l_ 
3a 



K 



247resceo^^ ' 
\{iB-2+ ^iB (35 - 4) 
7 96?™* 

^sc ~r 



2 



1/3 



/3F2 



(B16) 
(B17) 
(B18) 
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The density ndcpi can be calculated from the total band 
bending e^o = Ecb — Eqs (b=bulk, s=surface) , 



e?T-dopi = e {Na - No) z^^ ^2 {Na - N^) escfoe^i 



'Od, 



ecfiod = 600 - e(j)os - ksT. 



(B20) 
(B21) 



Here e^od is the band bending caused by the charges in 
the depletion layer (acceptors and donors) and e0os that 
caused by the free electrons in the inversion layer. To 
get this equations one has to assume that in the deple- 
tion layer all acceptors are charged and that in the bulk 
there is charge neutrality. The boundary between the de- 
pletion layer and the bulk is not sharp, this is described 
by the term —ksT in e^od-'^^ To calculate (f)Qs one has to 
solve the Poisson equation with the charge density which 
corresponds with the test wave function ip, the result is 



Sens 
besc^o ' 



(B22) 



As can be seen from Fig. 2] for the total band bending 



600 = Ecb ~ Ecs, 

600 = Ecb - + M 

600 = -fJ-Cb + A^-Eo - 



- Eq 



En ~ E, 



Cs, 



(B23) 



holds. The chemical potential relative to the conduction 
band edge in the bulk fj,cb is determined by the charge 
neutrality and can be calculated by solving 



}Mjb_ 



Na 



— — ^Cb 

ksT 



(B24) 



numerically. Nc and Ny are the effective densities of 
state in the conduction band and in the valence band, £a 
is the acceptor ionization energy, qa the acceptor degen- 
eracy factor, and Eg the gap energy. According to Sze^i 
this quantities are given by 



27r;i2 J 
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Nv^g. 



f nidhSBkBT 



3/2 



(B25) 
(B26) 
(B27) 



' ^^^^^^ T + /3si' 
The values we used can be found in Table HI 

The ground state energy Eq itself is a (small) part of 
the the total band bending, this problem is solved by a 
fix point iteration using = a-s start value. 
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